Measurement of Angular Distributions of Drell-Yan Dimuons in p + d Interaction at 

800 GeV/c 
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We report a measurement of the angular distributions of Drell-Yan dimuons produced using an 
800 GeV/c proton beam on a deuterium target. The muon angular distributions in polar angle 
6 and azimuthal angle 4> have been measured over the kinematic range 4.5 < m^^ < 15 GeY/c^, 
< pt < 4 GeV/c, and < xf < 0.8. No significant cos20 dependence is found in these proton- 
induced Drell-Yan data, in contrast to the situation for pion-induced Drell-Yan. The data are 
compared with expectations from models which attribute the cos2(^ distribution to a QCD vacuum 
effect or to the presence of the transverse-momentum-dependent Boer-Mulders structure function 
hi . Constraints on the magnitude of the sea-quark hj^ structure functions are obtained. 

PACS numbers: 13.85.Qk, 14.20.Dh, 24.85. +p, 13.88.-|-e 



The Drell-Yan process in which a charged lepton 
pair is produced in a high-energy hadron-hadron interac- 
tion via the qq —^ process, has been a testing ground 
for perturbative QCD and a unique tool for probing par- 
ton distributions of hadrons. The Drell-Yan production 
cross sections can be well described by next-to-leading 
order QCD calculations 0. This provides a firm theo- 
retical framework for using the Drell-Yan process to de- 
termine the antiquark content of nucleons and nuclei , 
as well as the quark distributions of pions, kaons, and 
antiprotons 

Despite the success of perturbative QCD in describing 
the Drell-Yan cross sections, it remains a challenge to un- 
derstand the angular distributions of the Drell-Yan pro- 
cess. Assuming dominance of the singic-photon process, 
a general expression for the Drell-Yan angular distribu- 
tion is 

^ oc 1 + Acos^ 6 + /I sin 20 cos -f ^ sin^ cos 20, (1) 

where 9 and denote the polar and azimuthal angle, 
respectively, of the in the dilepton rest frame. In 
the "naive" Drell-Yan model, where the transverse mo- 



mentum of the quark is ignored and no gluon emission 
is considered, A = 1 and /x = i/ = are obtained. 
QCD effects and non-zero intrinsic transverse mo- 
mentum of the quarks can both lead to A ^ 1 and 
/i, ^ 0. However, A and ly should still satisfy the rela- 
tion 1 — A = 2i/ [sj. This so-called Lam- Tung relation, 
obtained as a consequence of the spin- 1/2 nature of the 
quarks, is analogous to the Callan-Gross relation in 
Deep-Inelastic Scattering. While QCD effects can signif- 
icantly modify the Callan-Gross relation, the Lam- Tung 
relation is predicted to be largely unaffected by QCD 
corrections 9]. 

The first measurement of the Drell-Yan angular dis- 
tribution was performed by the NAlO Collaboration for 
TT- +W at 140, 194, and 286 GeV/c, with the highest 
statistics at 194 GeV/c 0. The cos 20 angular depen- 
dences showed a sizable p, increasing with dimuon trans- 
verse momentum (pt) and reaching a value of « 0.3 at 
Pt = 2.5 GeV/c (see Fig. 1). The observed behavior of v 
could not be described by perturbative QCD calculations 
which predict much smaller values of |6j • The Fermilab 
E615 Collaboration subsequently performed a measure- 
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FIG. 1: Parameters \,iJi.,v and 2i/ — (1 — A) vs. pr in the 
CoUins-Soper frame. Solid circles are for E866 p + d at 800 
GeV/c, crosses are for NAIO vr" + W at 194 GeV/c, and 
diamonds are E615 tt" + W a,t 252 GeV/c. The error bars 
include the statistical uncertainties only. 



ment oin^ +W Drell-Yan production at 252 GeV/c with 
broad coverage in the decay angle 9 11]. The E615 re- 
sults showed that A deviates from 1 at large values of 
(the Bjorken-x of the incident pions), and both fi and v 
have large non-zero values. Furthermore, the E615 data 
showed that the Lam- Tung relation, 2i^ = 1 — A, is clearly 
violated. (See Fig. 1.) 

The NAIO and E615 results on the Drell-Yan angu- 
lar distributions strongly suggest that new effects be- 
yond conventional perturbative QCD are present. Sev- 
eral attempts have been made to interprete these data. 
Brandenburg, Nachtmann and Mirke suggested that a 
factorization-breaking QCD vacuum may lead to a corre- 
lation between the transverse spin of the antiquark in the 
pion and that of the quark in the nucleon This would 
result in a non-zero cos 20 angular dependence consistent 
with the data. As pointed out by Boer et at, a possible 
source for a factorization-breaking QCD vacuum is helic- 
ity flip in the instanton model [13i] . Several authors have 
also considered hi ghe r-twist effects from quark- antiquark 
binding in pions jl4L llll . motivated by earlier work of 
Berger and Brodsky [ig. This model predicts behavior 



of n and v in qualitative agreement with the data. How- 
ever, the model is strictly applicable only in the — > 1 
region while the NAIO and E615 data exhibit nonpertur- 
bative effects over a much broader kinematic region. 

More recently, Boer pointed out (13 that the cos 2(f> an- 
gular dependences observed in NAIO and E615 could be 
due to the fc-r-dependent parton distribution function . 
This so-called Boer-Mulders function is an exam- 
ple of a novel type of fcr-dependent parton distribution 
function, and it characterizes the correlation of a quark's 
transverse spin and its transverse momentum, kx, in an 
unpolarized nucleon. It has an interesting property of be- 
ing a time-reversal odd object and owes its existence to 
the presence of initial/final state interactions fli|. The 
Boer-Mulders function is the analog of the Collins frag- 
mentation function j2(j|. which describes the correlation 
between the transverse spin of a quark and the trans- 
verse momentum of the particle into which it hadronizes. 
Model calculations for the nucleo n (p ion) Boer-Mulders 
functions have been carried out |2ll E3 | in the 

framework of quark-diquark (quark-spectator-antiquark) 
model, and can successfully describe the v behavior ob- 
served in NAIO (2^1 . 

To shed additional light on the origins of the NAIO and 
E615 Drell-Yan angular distributions, we have analyzed 
p + d Drell-Yan angular distribution data at 800 GeV/c 
from Fermilab E866. There are several physics motiva- 
tions for this study. First, there has been no report on the 
azimuthal angular distributions for proton-induced Drell- 
Yan - all measurements so far have been for polar angular 
distributions 0, Second, proton-induced Drell-Yan 
data provide a stringent test of theoretical models. For 
example, the cos2</) dependence is expected to be much 
reduced in proton-induced Drell-Yan if the underlying 
mechanism involves the Boer-Mulders functions. This 
is due to the expectation that the Boer-Mulders func- 
tions are small for the sea-quarks. However, if the QCD 
vacuum effect is the origin of the cos 2ip angular de- 
pendence, then the azimuthal behavior of proton-induced 
Drell-Yan should be similar to that of pion-induced Drell- 
Yan. Third, the validity of the Lam- Tung relation has 
never been tested for proton-induced Drell-Yan, and the 
present study provides a first test. 

The Fermilab E866 experiment was performed using 
the upgraded Meson-East magnetic pair spectrometer. 
Details of the experimental setup have been described 
elsewhere [2^. An 800 GeV/c primary proton beam with 
up to 2 X 10""^^ protons per beam spill was incident upon 
one of three identical 50.8 cm long cylindrical stainless 
steel target flasks containing either liquid hydrogen, liq- 
uid deuterium or vacuum. A copper beam dump located 
inside the second dipole magnet (SM12) absorbed pro- 
tons that passed through the target. Downstream of the 
beam dump was an absorber wall that completely filled 
the aperture of the magnet. This absorber wall removed 
hadrons produced in the target and the beam dump. 
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TABLE L Mean values of the X, ^1,1^ parameters and the quan- 
tity 22/ — (1 — A) for three Drell-Yan measurements. The pr 
dependence of these quantities is shown in Fig. 1. 





p + d 
800 GeV/c 
(E866) 


TV- +W 

194 GeV/c 
(NAIO) 


252 GeV/c 
(E615) 


(A) 


1.07 ±0.07 


0.83 ± 0.04 


1.17 ±0.06 




0.003 ± 0.013 


0.008 ±0.010 


0.09 ±0.02 




0.027 ±0.010 


0.091 ±0.009 


0.169 ±0.019 


(2z.-(l-A)) 


0.12 ±0.07 


0.01 ±0.04 


0.51 ±0.07 



Several settings of the currents in the three dipole mag- 
nets (SMO, SM12, SMS) were used in order to optimize 
acceptance for different dimuon mass regions. Data col- 
lected with the "low mass" and "high mass" settings 
on liquid deuterium and empty targets were used in this 
analysis. The detector system consisted of four track- 
ing stations and a momentum analyzing magnet (SMS). 
Tracks reconstructed by the drift chambers were extrapo- 
lated to the target using the momentum determined from 
the bend angle in SMS. The target position was used to 
refine the parameters of each muon track. 

From the momenta of the /i+ and fi^ , kinematic vari- 
ables of the dimuons {xp ,m^f^^pT) were readily recon- 
structed. The muon angles 9 and (p in the CoUins-Soper 
frame [ilj were also calculated. To remove the quarko- 
nium background, only events with 4.5 < m^^ < 9 
GeV/c^ or > 10.7 GeV/c^ were analyzed. A total 
of 118,000 p + d Drell-Yan events covering the decay an- 
gular range —0.5 < cos0 < 0.5 and —it < (jj < tt remain. 
Detailed Monte-Carlo simulations for the experiment us- 
ing the MRST98 parton distribution functions [l^l for 
NLO Drell-Yan cross sections have shown good agree- 
ments with the data for a variety of measured quantities. 

Figure 1 shows the angular distribution parameters 
A, /i, and v vs. pt- To extract these parameters, the 
Drell-Yan data were grouped into 5 bins in cos0 and 8 
bins in (p for each pT bin. A least-squares fit to the data 
using Eq. 1 to describe the angular distribution was per- 
formed. Only statistical errors are shown in Fig. 1. The 
primary contributions to the systematic errors are the 
uncertainties of the incident beam angles on target. The 
analysis has been performed allowing the beam angles to 
vary within their ranges of uncertainty. From this study, 
we found that the systematic errors are comparable to the 
statistical errors for each individual px bin. However, the 
Pt averaged values (A), (/i), and (v), are dominatd by the 
statistical errors. 

For comparison with the p + d Drell-Yan data, the 
NAIO TT- +W data at 194 GeV/c and the E615 Tr~ +W 
data at 252 GeV/c are also shown in Fig. 1. To test the 
validity of the Lam- Tung relation, also shown in Fig. 1 
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FIG. 2: Parameter v vs. pt in the Collins-Soper frame for 
three Drell-Yan measurements. Fits to the data using Eq. 3 
and Mc = 2.4 GeV/c^ are also shown. 



is the quantity, 2;^ — (1 — A), for all three experiments. 
For p + d a.t 800 GeV/c, Fig. 1 shows that A is consis- 
tent with 1 , in agreement with previous studies '3, ^23\ , 
while fi and v deviate only slightly from zero. This is in 
contrast to the pion- induced Drell-Yan results, in which 
much larger values of v are found. Table I lists the mean 
values of A, /i, v and 2i/ — (1 — A) for these three experi- 
ments. Again, the qualitatively different behavior of the 
azimuthal angular distributions for p + d versus 7r~ + W 
is evident. It is also interesting to note that while E615 
clearly establishes the violation of the Lam- Tung rela- 
tion, the NAIO and the p + d data are largely consistent 
with the Lam- Tung relation. 

In an attempt to extract information on the ma gnit ude 
of the hj^ function from the NAIO data, Boer [l2| as- 
sumed that hj^ is proportional to the spin-averaged par- 
ton distribution function /i: 



hi{x, k^) = Ch 



ar Mc'Mh 



h{^). (2) 



where fcr is the quark transverse momentum, Mu is the 
mass of the hadron H (pion or nucleon), and Mc and 
Ch are constant fitting parameters. A Gaussian trans- 
verse momentum dependence of e""^'^^ with ar = 1 
(GeV/c)~^ was assumed. The cos 20 dependence then 
results from the convolution of the pion hj^ / fi term with 
the nucleon / fi term, and the parameter 1/ is given as 



16ki 



p^^Ml 



- 4M2 )2 ' 



(3) 



where ki = ChiCh2/'^, and Hi, H2 denote the two inter- 
acting hadrons. As shown in Fig. 2, a good description 
of the NAIO data is obtained with ki = 0.47 ± 0.14 and 
Mc = 2.4 ± 0.5 GeV/c^. A fit to the E615 ly data at 
252 GeV/c using Mc = 2.4 GeV/c^, also shown in Fig. 
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FIG. 3: Parameter v vs. m^^, xp, xi, and X2 in the CoUins- 
Soper frame for p+d at 800 GeV/c. The error bars correspond 
to the statistical uncertainties only. 



2, gives Ki — 0.93 ± 0.10. These large values of ki sug- 
gest sizable hj^ functions for the valence antiquarks in 
the pion and for the valence quarks in the nucleon. 

A fit to the E866 p + d data using Eq. 3 yields ki = 
0.11 ± 0.04 for Mc = 2.4 GeV/c^, as shown in Fig. 2. 
As noted earlier, proton-induced Drell-Yan involves a va- 
lence quark annihilating with a sea quark. A comparison 
of the values of ki from proton-induced Drell-Yan with 
those from pion-induced Drell-Yan suggests that the ra- 
tio hj- 1 fi for the nucleon sea quarks is substantially be- 
low that for valence quarks. More specifically, the value 
of Ch for the sea is approximately a factor 4-8 smaller 
than that for valence quarks. 

The Drell-Yan angular distributions have also been an- 
alyzed for other kinematic variables. Figure 3 shows the 
values of v for p+d vs. runn, xf,xi, and X2, where xi and 
X2 are the Bjorken-x for the beam and target partons, 
respectively. Again, for each bin the data were divided 
into 5 bins in cos^ and 8 bins in <p in order to extract 
the angular distribution parameters. Figure 3 shows no 
significant dependence on these kinematic variables. 

In summary, we report a measurement of the angu- 
lar distributions of Drell-Yan dimuons for p + at 800 
GeV/c. The pronounced cos2</) azimuthal angular de- 
pendence observed previously in pion-induced Drell-Yan 
is not observed in the p + d reaction. The Lam- Tung 
relation, found to be strongly violated in the E615 pion- 
induced Drell-Yan, remains largely valid for p + d Drell- 
Yan. These results put constraints on theoretical mod- 
els that predict large cos 2<p dependence originating from 



QCD vacuum effects. They also suggest that the Boer- 
Mulders functions hj- for sea quarks are significantly 
smaller than those for valence quarks. 

This work was supported in part by the U.S. Depart- 
ment of Energy and the National Science Foundation. 
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